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Abstract: We present a general method for obtaining a lower bound for the ground state 
entropy density of the Ising Model with nearest neighbor interactions. Then, using this 
method, and with a random coupling constant configuration, we obtain a lower bound for 
the ground state entropy density of the square, triangular, and hexagonal two-dimensional 
lattices with free, cylindrical, and toroidal boundary conditions. 



1 Introduction 



1.1 Overview 

h ; 

In this paper, we establish a general method for obtaining a lower bound for the ground 
state entropy density So of arbitrary lattices with arbitrary boundary conditions. We then 
use this method on two-dimensional square, triangular, and hexagonal lattices with free, 
cylindrical, and toroidal boundary conditions — while assuming a random coupling constant 
configuration — to obtain numerical values for these lower bounds. Let us now note that 
our method to be presented is derived from a method employed by Avron et al [1], in 
which a value of approximately 1.100 x 10~ 19 is obtained as a lower bound for Sq for the 
two-dimensional square lattice with a random coupling constant configuration. Using our 
method, we obtain a value of approximately 4.883 x 10 -6 for this same lattice. Let us now also 
note that Loebl and Vondrak [2] obtained an upper bound of 0.548 for So for a certain class 
of square lattices with toroidal boundary conditions. Loebl and Vondrak then conjecture 
that this upper bound holds for all square lattices with toroidal boundary conditions. If this 
is indeed true, then coupled with our result, we find that the ground state entropy density 
of the square lattice with toroidal boundary conditions and a random coupling constant 
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configuration is bounded by 



4.883 x KT 6 < Sn < 0.548. 



However, it should be noted that while our lower bound of 4.883 x 10~ 6 is many orders 
of magnitude greater than the previous lower bound of 1.100 x 10~ 19 , it is still far from 
what one would expect based on the upper bound of 0.548. This discrepancy is most likely 
attributable to the method which we use in this paper to establish the lower bound. 

1.2 Background Information 

The Ising Model is defined by a finite set A, representing atoms, and a coupling constant 
configuration J : B(A) — > {—1,1}, describing the interaction between neighboring atoms. 
The set B(A) is the set of all bonds in A. Formally, we write B(A) C {{si,Sj} : Si,Sj G 
A,i 7^ j}. We call elements of A sites, and often denote the bond between the sites Si,Sj 
by bij. Note that with this definition of the Ising Model, we are not restricting ourselves to 
only subsets of Z 2 , but in fact are considering arbitrary graphs. A spin state for A is an 
assignment of either spin up or spin down to each atom; that is, a spin state is a function 
a : A — > { — 1,1}. The energy of the system is given by 



The ground state energy is the minimum possible energy. Let us now note that if Jija^aj = — 1 
for some bond b^, then this in fact increases the energy; hence, we call such a bond unhappy. 
The set of all unhappy bonds in A with J and a is denoted U(A, J, a). Given A with a J, 
a spin state a that yields the ground state energy is called a ground state. The set of all 
ground states is denoted by -Do (A, J), and the number of ground states is called the ground 
state degeneracy. We will henceforth be considering only connected sets A with a given J. 
However, to properly define connectedness, we must first introduce the notion of a curve. 

Definition 1.1 Given A, a curve of length n > 1 in A is a finite sequence s : {0, 1, . . . , n} — > 

S, where S C A and {sj_i,Sj} G B(A) for all 1 < i < n. The curve s is said to connect 
the sites s an d s n , and the set of all curves connecting s and s n is denoted by C(so,s n ). 
Finally, we call s a closed curve iff s n = s . 

The notion of a closed curve will become useful in the next section. We can now say 
what me mean by connectedness. By a set A being connected, we mean that for any two 
sites x, y G A, there exists a curve c G C(x,y) connecting them. 

In this paper, we obtain a lower bound for the ground state entropy density; to do this, 
however, we must first obtain a lower bound for the ground state degeneracy, which in turn, 
requires the following three preliminary results to establish our formalism. 

2 Preliminary Results 

Our first result is the following trivial theorem which tells us that in order to minimize the 
energy we need to minimize the number of unhappy bonds. 




bij€B(A) 
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Theorem 2.1 Given a spin state a for A, H(A, J, a) = 2\U(A, J, a)\ - |S(A)|. 
Proof: Letting U = U(A, J, a), we have 

H(A,J,a) = — 2J JijO r i&j = —\ 2J JijCTi(Tj+ 2J Jij&i&j 

bij€B(A) \6y£B(A)nf7 b i:j eB(A)\U 

E -1+ E 1 

v 6ijeB(A)nl/ b i:j €B(A)\U 

-(-\U(A,J,a)\ + \B(A)\-\U(A,J,a)\) 
2\U(A,J,a)\-\B(A)\. 



The next observation is a simple one, but fundamental to understanding our formalism. 

Theorem 2.2 Given a spin state a for A. Then, for any closed curve s in A of any length 
n, we have 



ra—l n— 1 



i=0 i=0 

Proof: This is true since 

n-1 



j J Ji,i+\0~iO~i+l = (</oi (T (T l)(^12 cr l cr 2) ' • • • ' ( J„-2,n-l cr n-20"n-l) {Jn-l,n°~n-l°~n) 
i=0 

fn-l \ /n-1 

n^ + i n 



a 2 



,i=0 / \i=0 

n-1 



i=0 



Theorem |2.2| tells us how the distribution of unhappy bonds depends on J: along each 
closed curve in A, the parity of the number of negative J bonds equals the parity of the 
number of unhappy bonds. This observation will be used extensively in our method of 
bounding the ground state degeneracy. 

We are now ready for our third preliminary observation which treats the following general 
question. Given a spin state a for A, what are all the spin states that give the lattice the 
same energy as cr? 

To begin, let S be any subset of A, and let erg be the spin state for A obtained from a by 
flipping the spins of the sites in S. How do the lattice energies H (A, J, a) and H (A, J, as) 
compare? 

This question can be answered by defining the set Bs = {hj G B(A) : s, G S, Sj ^ S}. 
The set Bs can be easily visualized by a contour, as in Picture 1, where we have only shown 
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the special case that A is a subset of a square lattice. Note that by flipping the spins of only 
the sites in S, only the happiness of bonds in the set B s change. Indeed, for if the spins 
of Si and Sj are both flipped, the happiness remains unchanged. Moreover, if neither the 
spin of S{ nor Sj is flipped, then the happiness trivially remains unchanged. Thus, letting 
U = U(A, J, a), to get the set Us = U(A, J, as) of unhappy bonds in A with J and as, we 
simply "complement" the bonds along the contour; that is, 

U s =[U\(B s nU)]U(Bs\U). 

From this, it immediately follows that \U$\ = \U\ — \B S HU\ + \B S \U\, which coupled with 
Theorem |2.1| tells us how the lattice energies compare. We have thus proved the following 
theorem. 

Theorem 2.3 Given A with J and a, let U = U(A, J, a). For any S C A, define as and 
Us as above. 

Then, H{A, J, a) = H(A, J, a s ) iff\B s nU\ = \B S \ U\ . 

Picture 1 
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The contour is shown by the thick line. The sites of S 
are those enclosed by the contour. The set Bs is the set 
of bonds that cut the contour perpendicularly. 

Thus, to find the set of all spin states that give the lattice the same energy as a, we need 
to find all subsets S of A such that the number of bonds in Bs that are in U is equal to the 
number of bonds in Bs that are not in U; that is, \B$ H U\ — \B$ \ U\. We call such a set S 
that satisfies \B$ Pi U\ — \Bs \ U\ entropic relative to U, or entropic relative to a. 

3 Main Results 
3.1 Local Degeneracy 

We can now discuss the idea behind the paper. We construct a J configuration in a subset 
M of A that guarantees a local degeneracy. By this we mean that we are guaranteed to find 
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degeneracy within this subset independent of the J configuration outside of this subset. We 
call such an entity a module. We then suppose to have a random J configuration for A, and 
then using probabilistic methods we determine the number of expected modules in all of A, 
thereby obtaining a lower bound for the ground state degeneracy. 
We now proceed to rigorously define what a module is. 

Definition 3.1 A module in A is a set MCA with a Jm such that for any J for A and 

any ground state a G D (A,J), if Jm = J\b(m), then there exists a nonempty set S C M 
such that a s G D (A, J). 

That is, we take a set of sites M in A and define what the coupling constant configuration 
for it is. We then require M to have the property that given any ground state a, there exists 
at least one non-trivial entropic set relative to a contained in M. From this it immediately 
follows that if we have a collection of n pairwise disjoint modules in A, then a lower bound 
for the ground state degeneracy for A with any J is \D (A, J)| > 2 n . 

Let us now suppose the existence of a module M. We will now demonstrate the method 
by which we obtain a lower bound for the ground state degeneracy. We first let k be the 
maximum number of pairwise-disjoint subsets T of A, where each Tj has the possibility of 
being the module M. Note that if {Tj}^ =1 is a partition of A, then k = j^.. We now let p be 
the probability that an arbitrary bond is a negative J bond. Then, the probability of any 
Tj being the module M will be some function of p that depends on M. We will denote this 
probability by f(p)- We also want to consider site and bond diluted lattices. To understand 
diluted lattices, it is helpful to begin with a "first" lattice A x with the set of bonds B(A 1 ), 
and two probability functions P s and P b — P s defined on Ai and P b defined on -B(Ai). With 
this set up, a site diluted lattice is any subset A of Ai that is chosen consistent with P s , 
where the interpretation of P s is that P s (x) = P[x G A]. A site diluted lattice can be used to 
model a substance that has sporadic impurities, but that is homogeneous otherwise. Semi- 
conducting silicon wafers doped with boron is one example. Next, A is a bond diluted lattice 
if its set of bonds B(A) C B(A 1 ) is chosen consistent with P b , where the interpretation of 
P b is that P b (x) = P[x G B(A)]. A bond diluted lattice can be used to model an imperfect 
crystal; here, there are some sets {s, s'} of neighboring atoms that interact in the perfect 
crystal Ai, but that do not interact in the imperfect crystal A, which is in fact what makes 
A imperfect. We will suppose that P s and P b are constant functions, and denote their values 
by p s and p b , respectively. With these definitions, it now follows that the probability that an 
arbitrary Tj is the module M is given by f(p)p[ M ^p b B< " M ^ ■ With this in mind, we now prove 
the first main result of our paper. 

Theorem 3.2 Given an infinite set of sites T with a J, suppose that {T i }^ 1 is a pairwise- 
disjoint collection of subsets of T such that each % has the possibility of being the module 
M . Then, for any e > and for any 5 > 0, there exists a ko such that for all k > ho, if 

k 

ACT and Q Tj C A, then 



i=l 



\Dq(A, J\b(a)) I > 2 



*/(p)pl M| pi S(M)l (l-e) 



> 1 
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Proof: Let e > and 5 > be given. Since the probability of any Tj being the module M 
is /(p)pL M 'f>!, B ^', the Law of Large Numbers furnishes a k G N such that for all k > k , 



where n is the number of the first k Tj's that are modules in T. Now let A C T be given 

k 

such that |J Ti C A. Then, there are at least n pairwise disjoint modules M in A, and hence 

i=l 

| -D (A, J\b(A))\ > 2 n . The conclusion now follows. a 
3.2 Existence of Modules 

We next wish to demonstrate the existence of modules for square, triangular, and hexagonal 
lattices. Then, using Theorem pT2| , we will calculate explicitly a lower bound for each of these 
lattice types with the restriction that we have a random J configuration — the probability 
that a bond is a negative J bond is \ — and that the lattice is neither bond nor site diluted. 
That is, p s — Pb — 1- But we will first need to calculate /(|) for each of the modules. So 
as to make the calculation of /(|) a trivial task, we must introduce the notion of plaquettes 
and frustrated plaquettes. If we have a lattice A, a plaquette is the basic building block of 
this lattice. For example, for a square, triangular, or hexagonal lattice, a plaquette is a set of 
four, three, or six sites, respectively, that are the vertices of the smallest square, triangle, or 
hexagon, respectively. Note that a plaquette can be realized as the range of a closed curve. 
Now suppose we have a plaquette in which the product of the J's along this closed curve is 
— 1. Then, by Theorem [2.2|, no matter what spin state is put on A, there must be at least 
one unhappy bond in this plaquette. We therefore call such a plaquette in which the product 
of the J's is —1 a frustrated plaquette. With this terminology and the following lemma, we 
will be able to compute f(p) for random J's. 

Lemma 3.3 Let s be a closed curve of length n with q < n of its J values specified. Then, 
the number of ways to get an odd number of negative J bonds along s is equal to the number 
of ways to get an even number of negative J bonds along s. 

Proof: Without loss of generality, suppose that of the q bonds already specified, an odd 
number are negative J bonds. We can also without loss of generality assume that n — q is 
even. Then, the number of ways to get an odd number of negative J bonds is 





The number of ways to get an even number of negative J bonds is 




However, it is a well known fact that these two are equal. 
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Given a random J configuration, this lemma shows us that as long as there is at least one 
J value not specified on a plaquette, then the probability that this plaquette is frustrated 
is |, and the probability that it is unfrustrated is ~. Thus, given a set of plaquettes whose 
frustrations are specified on exactly m of these plaquettes, it follows that the probability 
of this frustrated plaquette configuration occurring is (|)"\ provided of course that we can 
snake around in such a way that we have at least one degree of freedom for each plaquette. 
Therefore, /(|) = (§) m - 

Having calculated /(|), we now proceed to prove the existence of modules. We begin by 
proving the existence of a module for a square lattice. 

Theorem 3.4 Let A be a square lattice and suppose M C A is the set of '25 sites with any 
Jm that gives the plaquettes Pi, ■ ■ ■ ,pu the frustrations as indicated in Picture 2. Then, M 
is a module in A. 
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In pictures 2-4, numbers at vertices refer to site numbers and num- 
bers in the centers of plaquettes refer to plaquette numbers. If a 
plaquette number is circled, then that plaquette is specified as frus- 
trated; otherwise it is specified as unfrustrated. The frustration 
of all plaquettes that do not have a number is left unspecified. 



Proof: Suppose a Q is a ground state for A with any J such that J\b(m) — J Mi and suppose 
that M is not a module in A. Then, if S C M is nonempty, it follows that 5* is not entropic 
relative to U = U(A, J, a ). 

Define M' = {si, . . . , S23}. By Theorem [E3] we know that any site Sj G M' can be 
involved in at most one unhappy bond in Uq] for if Sj were involved in more, then the set 
{s^ would be entropic relative to U . Now consider the plaquette p 5 . Since p 5 is frustrated, 
we know that there is at least one bond in p 5 that is also in Uq. We thus have four cases: 

Case(l) b 6J G U 

By repeated application of Theorem |2]^ to the plaquettes p%, P2, P3, P7, pn, Pw, P9, in 
that order, we deduce that the following bonds are in Uq: 61,2, ^3,8, ^4,9, &13,H, &18,19> ^12,17, 
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611,16, as seen m Picture 2. We now observe that the set {sj, s 8 , Si 2 , Si 3 } is entropic relative 
to U which is a contradiction. Therefore, Case (1) is not possible. 

Following the structure of the reasoning in Case (1), we similarly find that none of the 
cases (2) 6 6j n 6 Uo, (3) &n ; i2 G Uq, or (4) 6712 6 Uq is possible. But this contradicts the fact 
that at least one of these cases must hold. Therefore, M is a module in A. B 

Before we continue, there are two important things we need to be clear about. First, 
thus far we have only shown that if there is a Jm that yields the configuration of frustrated 



plaquettes in Picture 2, then M is a module in A. To be able to apply Theorem 3.2, we 



need to further show the existence of such a Jm- As with each of our modules, however, the 
existence of such a Jm is a trivial fact, and we leave its verification to the reader. Second, 
to apply Theorem 3.2 with p = ^, we need to verify that for each of our modules, we can 



snake around the module in such a way that we have at least one degree of freedom for each 
plaquette in the module. This too is an easy exercise that we leave to the reader. 

Having said this, we next observe that if we rotate the module in Picture 2 by 90 degrees, 
then by symmetry, this will be another module. We thus have two modules for a square 
lattice. 

We now note that we can define a subset A of a square lattice that can be partitioned 
by 5 x 5 subsets, and in such a way that A can have either free, cylindrical, or toroidal 
boundary conditions. Thus, since both of our modules have 25 sites and we have specified 
the frustration of 14 plaquettes in each, it follows that the probability that any given 5x5 
subset in the partition is one of our modules is 2kf{\) = 2 ^ 13 . Thus, we can take a number 

|A| 



arbitrarily close to 2 5 ' as a lower bound for the ground state degeneracy of a square 
lattice. 

We now demonstrate the existence of a module for a triangular lattice. The module is 
seen in Picture 3. The proof that this is a module closely resembles the previous proof. 



Theorem 3.5 Let A be a triangular lattice and suppose M C A is the set of 21 sites with 
any Jm that gives the plaquettes Pi, ■ ■ ■ ,pig the frustrations as indicated in Picture 3. Then, 
M is a module in A. 

Proof: Let o"o be a ground state for A with any J such that J\b(m) — Jm, and suppose that 
M is not a module in A. Define M' = {s±, . . . , su}. Then, any site in M' can be involved in 
at most two unhappy bonds. Since plaquette pn is frustrated, it has at least one unhappy 
bond. By symmetry, we can assume this bond to be 6 56 . We then see that plaquette p 5 
needs one more unhappy bond since it is unfrustrated. By symmetry, we can take this bond 
to be 625- Next, since p^ is unfrustrated and ^3 is frustrated, we see that the bonds bu and 
614 must be unhappy. We then see that plaquette p\ is unfrustrated. Since it already has 
one unhappy bond, it must have one more. We therefore have an entropic set consisting of 
either site s\ or s 2 > which is a contradiction. Thus, there is indeed a non-trivial entropic set 
in M relative to ctq, and hence M is a module in A. B 

We now observe that we can define a subset A of a triangular lattice that can be parti- 
tioned by M, and in such a way that A can have any of the three boundary conditions we are 
considering. Thus, since \M\ = 21 and we have specified the frustration of 19 plaquettes in 
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Picture 3 




M, we see that for a random J, a lower bound for the ground state degeneracy of a triangular 

lattice can be taken arbitrarily close to 2 212 . 

Let us now turn our attention to hexagonal lattices. 

Theorem 3.6 Let A be a hexagonal lattice and suppose M C A is the set of 54 sites with 
any Jm that gives the plaquettes pi, . . . ,pig the frustrations as indicated in Picture 4- Then, 
M is a module in A. 

Proof: Let cr be a ground state for A with any J such that J\b(m) — Jm, and suppose 
that M is not a module in A. Letting Uq = U(A, J, er ), we first note that if S is any subset 
of M, then \B S \U \ > \B S C\Uo\ since a is a ground state. Also, since each site s G M is 
involved in at most three bonds, it follows that at most one of the bonds containing s is in 
U . 

Now consider plaquette p\. Since pi is frustrated, at least one of the bonds in p 1 must 
be in U . First consider the case that 6 2 3 £ U . 

We now observe that if &io,n G U , then ^io^i = {§2, Sio} is entropic relative to U , while 
if &i2,i3 G Uo, then 5'i2,i3 = {«3, S13} is entropic relative to Uq. In view of this, and the fact 
that P3 is unfrustrated, we conclude that 611,12 G Uq. 

Using reasoning similar to that just used, together with the fact that pi 2 is frustrated, 
we see that either b u , 37 G U or 6 36i37 G U . If 614,37 e ^0, then ^14,37 = {s 3 , s 12 , s 13 , s 14 ] is 
entropic relative to U . Hence, 636,37 £ ^o- 
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Picture 4 




From this and the fact that pi 3 is frustrated, we either have 615,40 G Uq or 6 3 g,4o G t/o- If 
615,40 G C/ , then 515,40 = {s 3 , S12, s 13 , s u , s 15 , s 37 } is entropic relative to U Q . Thus, 6 39 , 40 G C/ . 

Consequently, since pu is frustrated, either 616,17 G C/o or 617,42 G Uq. If 616,17 G C/o, then 
defining S W:17 = {s 3 , s 4 , si 2 , Si 3 , s i4 , S15, «ie, s 37 , s 40 } shows that |-B Sl617 \ C/ | < |-B Sl6il7 n U \, 
which is a contradiction. Thus, 617,42 G C/o- 

However, we now observe that 5*17,42 = Si 6 ,i7 U {S17} is entropic relative to Uq, and so we 
conclude that 617,42 ^ Uo, which is a contradiction. 

Therefore, &2 3 ^ Uo, and hence, by symmetry, none of the bonds in pi is in Uq, contra- 
dicting the fact that p\ is frustrated. 

Therefore, M is a module in A. B 

As with the triangular lattice, we see that we can define a subset A of a hexagonal lattice 
that can be partitioned by M and that can have any of the boundary conditions we are 
considering. Therefore, since \M\ = 54 and we have specified the frustration of 19 plaquettes 
in M, we see that for a random J, a lower bound for the ground state degeneracy for a 

hexagonal lattice can be taken arbitrarily close to 2 54 2 . 

4 Conclusion 

We now summarize our results of the lower bounds for the various lattices assuming a random 
J and no dilution of sites or bonds. Note that in every case, a lower bound is of the form 
2 c|A|(1 e) for some positive constant c. This constant c is in fact a lower bound for the ground 
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state entropy density of the lattice A, that is, the ground state entropy per unit volume, 
defined by 

g.(A,J) = '° fe(l ^ A - J)l) . 

We therefore have the following lower bounds for the ground state entropy densities of square, 
triangular, and hexagonal lattices, with a random J, and with either free, cylindrical, or 
toroidal boundary conditions: 
For a square lattice we have 

4.883 x KT 6 . 



25 • 2 13 

For a triangular lattice we have 

1 



c = 



21 • 2 19 

For a hexagonal lattice we have 

1 



c = 



54 • 2 19 



9.083 x 1(T 8 . 



3.532 x 1(T 8 . 



We should now note that the ground state entropy density S — as calculated via the 
Ising Model — is not necessarily the same as the ground state entropy density observed in the 
laboratory. Indeed, the temperature needed to put a system into its ground state may be 
unattainable in the laboratory. Hence, in the laboratory we take the thermodynamic limit 
first, and then we take the T — > limit. Let us denote this "observed" ground state entropy 
by S. However, the Ising Model is a model at zero temperature. Thus, S is obtained by 
taking the limit T — > and then the thermodynamic limit. This interchange of limits is 
cause for concern. However, Aizenman and Lieb [3] showed that 

S = lim sup So, 

A|oo 

where the supremum is taken over all boundary conditions. Therefore, the values of c we 
obtained for the various lattices are indeed lower bounds for the observed ground state 
entropies in the laboratory — the physically meaningful quantity. 

In closing, we emphasize that the formalism presented here is in no way restricted to the 
specific lattices or boundary conditions we have considered. Rather, it is a general method 
to find a lower bound for the ground state degeneracy of an arbitrary lattice with arbitrary 
boundary conditions. 

Also, while we didn't make explicit reference to the concept of a "dual lattice" or mini- 
mizing "path length" to minimize energy, it is these concepts that began our investigation, 
offered much insight into the problem, and led to the formalism given here. We feel that it 
would be very helpful to become acquainted with these concepts in applying this formalism, 
and for this purpose we refer the reader to [4] and [5], where these concepts make one of 
their earliest appearances in the literature. 
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